Bethe-Salpeter equation: 3D reductions, heavy mass limits 

and abnormal solutions. 



J. Bijtebier* 

Theoretische Natuurkunde, Vrije Universiteit Brussel, 
Pleinlaan 2, B1050 Brussel, Belgium. 
Email: jbijtebi@vub.ac.be 

February 9, 2008 



Abstract 

We show that the 3D reductions of the Bethe-Salpeter equation have the same bound state spectrum as the original 
equation, with the possible exception of some solutions for which the corresponding 3D wave function vanishes. 
The abnormal solutions of the Bethe-Salpeter equation (corresponding to excitations in the relative time-energy 
degree of freedom), when they exist, are recovered in the 3D reductions via a complicated dependence of the final 
potential on the total energy. We know however that the one-body (or one high mass) limit of some 3D reductions 
of the exact Bethe-Salpeter equation leads to a compact 3D equation (by a mutual cancellation of the ladder and 
crossed graph contributions), which does not exhibit this kind of dependence on the total energy anymore. We 
conclude that the exact Bethe-Salpeter equation has no abnormal solution at this limit, or has only solutions for 
which our 3D wave function vanishes. This is in contrast with the results of the ladder approximation, where no 
such cancellation occurs. We draw the same conclusions for the static model, which we obtain by letting the mass 
of the lighter particle go also to infinity. These results support Wick's conjecture that the abnormal solutions are 
a spurious consequence of the ladder approximation. 

PACS ll.10.Qr Relativistic wave equations. 

PACS 11. 10. St Bound and unstable states; Bethe-Salpeter equations. 
PACS 12.20.Ds Specific calculations and limits of quantum electrodynamics. 

Keywords: Bethe-Salpeter equations. Salpeter's equation. Breit's equation. 
Relativistic bound states. Relativistic wave equations. 



1 Introduction. 

A well known difficulty of the Bethe-Salpeter equation |], ||] is the existence (for some interactions at least) of abnormal 
solutions, corresponding to excitations in the relative time-energy degree of freedom || [|, |J ||. Wick suggested 
that these solutions could be an artefact due to the ladder approximation, and would disappear if the higher-order 
irreducible graphs could be included || ■ Mugibayashi showed however that the abnormal solutions remain present in 
the exact static model , which is simple enough to be completely solved. Recently, the conclusions of Mugibayashi 
were questioned by Jallouli and Sazdjian ||, so that the road to Wick's interpretation is open again. This possible 
conspiration of the higher-order contributions for suppressing the abnormal solutions reminds us a similar cancellation 
of the higher-order terms at the one-body (or one large mass) limit of some 3D reductions of the Bethe-Salpeter 
equation, leaving only a simple Dirac or Klein-Gordon equation |j, |l0|, [ll], |l2|, 13, A way of relating these two 



facts can be found in our previous study of a zero-D reduction of the ID Bethe-Salpeter equation for the static model 
at the ladder approximation (T^j. We showed that the reduction of the Bethe-Salpeter equation does not in general kill 
the possible abnormal solutions. Despite the elimination of the relative time, they are recovered via the multiplicity of 
solutions due to the complicated dependence of the final potential on the total energy. This was the starting point of 
the present work. It can be easily verified that this dependence of the reduced equation potential on the total energy 
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is a necessary condition to the presence of abnormal solutions of any reduction of a Bethe-Salpeter equation. As this 
dependence disappears from some well chosen 3D reductions of the exact Bethe-Salpeter equation at the one-body 
limit (and also from the zero-D reduction of the exact static model obtained by letting the mass of the lighter particle 
go also to infinity), the 3D equations have no abnormal solution at this limit. Conversely, these abnormal solutions 
are present in the reduction of the static model Bethe-Salpeter equation at the ladder approximation, and we do not 
see why they would not be a fortiori present in the case of finite masses. 

Does the absence of abnormal solutions for a 3D reduction imply their absence for the original Bethe-Salpeter 
equation itself? We shall show that the existence of a solution of the Bethe-Salpeter equation for a given total 
energy implies in general the existence of a solution of the reduced equation for this same total energy (there could 
be convergence problems, in principle manageable by analytic extensions). However, the reduction being based on 
a constraint fixing the value of the relative energy, a solution of the Bethe-Salpeter equation will not lead to a 
corresponding solution of the reduced equation if it happens to vanish for this value. This accident does clearly not 
happen with the normal solutions of the exact Bethe-Salpeter equation, and, from the study of the static model, we 
can deduce that it does not happen with at least some of the abnormal solutions of the ladder approximated Bcthc- 
Salpeter equation. For the possible abnormal solutions of the exact Bethe-Salpeter equation, we have of course no 
example to examine (as we are just trying to demonstrate that such solutions do not exist!). We can only conclude 
that these solutions, if they exist and survive at the one-body limit, are not transmitted to the 3D reductions we use. 

In section 2, we present our 3D reduction of the exact or ladder approximated Bethe-Salpeter equation for two 
fermions, the existence condition for this reduction and how to invert it. In section 3, we compute the one-body 
limit of this reduction in the QED case, the one-body limit of the Bethe-Salpeter equation itself and verify that the 
reduction of the limit is the limit of the reduction. In section 4, we get the static model by letting the mass of the 
former light fermion become also infinite, and compare again the zero-D reduction of the zero-body limit with the 
limit of the reduction. We also verify that the existence condition for the reduction is satisfied at least by the solutions 
of the ladder approximated Bethe-Salpeter equation for the static model. In section 5, we consider the Bethe-Salpeter 
equation and its 3D reduction for equal masses, and get the static model directly by letting this common mass go to 
infinity. This leads to the same Bethe-Salpeter equation as in the preceding section, but to another zero-D reduction, 
based on a different choice of the fixed relative energy. Section 6 is devoted to conclusions. 



2 The Bethe-Salpeter equation for two fermions and its 3D reduction. 

We shall write the Bethe-Salpeter equation for the bound states of two fermions as: 

$ = G A'$, (1) 

where $ is the Bethe-Salpeter amplitude, function of the positions x\,x 2 or of the momenta pi,p 2 of the fermions, 
according to the representation chosen. The operator K is the Bethe-Salpeter kernel, given in a non-local momentum 
representation by the sum of the irreducible two-fermion Feynman graphs. The operator Go is the free propagator, 
given by the product G01G02 of the two individual fermion propagators: 



Goi — — . - Pi — — 5 2 , • Pi (2) 

PiO — hi + ieh t -m\ + le 



where the hi are the Dirac free hamiltonians 

hi = di .fi + f3 t m l (i = 1, 2). (3) 

The Bethe-Salpeter kernel A should contain charge renormalization and vacuum polarization graphs, while the 
propagators Goi should contain self-energy terms (which can be transferred to K 16, In this work, we consider 
only the free fermion propagators in Goi and the "skeleton" graphs in A, and we hope that the inclusion of the 
various corrections would not change our conclusions. We shall define the total (or external, CM, global) and relative 
(or internal) variables: 

X = -(x l +x 2 ), P = Pi+p 2 , (4) 

x = xi-x 2 , P=-z(pi-P2), (5) 
and give a name to the corresponding combinations of the free hamiltonians: 

S = h 1 + h 2 , s = - h 2 ). (6) 
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In order to keep our formula's compact, we do not write the arguments of the functions and operators. When needed, 
we work in the momentum representation (most often, we indicate only the relative energy explicitly). 



2.1 3D reduction with Breit's propagator. 

The free propagator Go will be approached by a carefully chosen expression Gs, combining a constraint fixing the 
relative energy, and a global 3D propagator. The argument of the constraint and the inverse of the propagator 
should be combinations of the operators used in the free equations (at last approximately and for the positive-energy 
solutions). We shall choose Breit's propagator jl8j 

G 5 = -2in - J— 5(p - n)(3 x (3 2 (7) 

with 

M= — (Ei -El) = — (8) 

^ = ^=(p1 + m 2 )*. (9) 
We shall write the free propagator as the 3D approached propagator, plus a rest: 

G Q = Gs + G R . (10) 
The Bethe-Salpeter equation becomes then the inhomogeneous equation 

$ = G K<S> = (G s + G R )K<5> = * + G R K $, (11) 

with 

* = G S K<S> (= GsGq 1 ®). (12) 
Solving (formally) the inhomogeneous equation ( |TT| ) and putting the result into (|l2|), we get 

* = G s K(l - G R K)-H = G S K T ^ (13) 

where 

K T = K(l - GrK)- 1 = K + KG r K + ... = (1 - KGr^K (14) 

obeys 

K T = K + KGrKt = K + K T G R K. (15) 

The reduction serie s fll4| ) re- introduces in fact the reducible graphs into the Bethe-Salpeter kernel, but with Go replaced 
by Gr. Equation ((13|) is a 3D equivalent of the Bethe-Salpeter equation. It depends on the choice of Gs and is not 
manifestly covariant (of course, we can always write it in the center of mass frame (assuming P 2 > 0) and then replace 
everything by its covariant equivalent). We can write it as a pair of coupled equations by multiplying it by Pq — S 
and po — (Lt: 

(P - S)V = -2m8{p - fi)^ 2 K T ^, (16) 
(p - tx)V - 0. (17) 
The relative time can be completely eliminated by making 

* = V2^S(p a -ti)^' (18) 

which gives 

(Po - 5)*' = V'V' (19) 

V' = -2^/3^2 K T (fx, (j,). (20) 

The Po spectrum of (|l9| ) should be that of the original Bethe-Salpeter equation (|l|), without the possible states 
for which the expression used in (|l^) to define the wave function ^ vanishes. The normal bound state spectrum, 
connected to the nonrelativistic spectrum, will be obtained at first order by keeping only the first term K, (or an 
approximation of it) in Kt, the remaining of Kt being treated as a perturbation. If the Bethe-Salpeter equation 
has also an abnormal spectrum, it should, at least implicitly, remain present also in the 3D reduction (|l9|), via the 
Po— dependence of V . Writing 

Po*' = [S + V'(Pb)]V (21) 
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and diagonalizing the hamiltonian 

[S + V'(P )] %(P ) = ^(Po)^-(Po) (22) 
we get the energy spectrum as the set of the roots of the equations 

P -Ei{P ) = 0. (23) 

For each value of i, we should find the energy of the normal solution, plus possibly a set of abnormal energies. This 
scheme is in general too complicated to be used in actual numerical calculations. We know three exceptions: 
— When we want only to compute the normal solutions by perturbations. 

In the static model for two fermions in QED at the ladder approximation (one photon exchange) jl5). In this 
case, (S + V) is not an operator anymore, it is simply a function of the energy and the coupling constant. The 
reduced equation (a simple numerical equation) gives the only normal solution and the even abnormal solutions (for 
the odd abnormal solutions, ^ vanishes, as the relative energy is fixed to zero in the reduction performed in jl5| ). 
The unique function Pi(Po) is divided into sectors separated by poles and containing each one solution. The reduction 
series converges only in the first sector containing the normal solution, but can be analytically extended by using the 
Pade approximants method. 

— At the one-body or zero-body limits of the exact equation, since the final 3D potential reduces then to a simple 
compact term, as we shall see in section 3. 

2.2 Existence conditions and inversion of the reduction. 

The 3D reduction performed above becomes impossible when the 3D wave function 'J, defined by ((l^), is the null 
function. This gives the condition 

* = G S K<$> ^ -> (K$)(p) (24) 
or, using the Bethe-Salpeter equation itself: 

*' = TT~q G o Hm) $ (M) = ^#( J Po-5)(P 2 -4 s 2 )$( M ) jt 0. (25) 

The expression 

(P - S)(P 2 - 4s 2 ) = (P -ht- h 2 )(P -hi + h 2 ){P Q + tn- h 2 ), (26) 

when projected on the subspaces hi — ±-Ej, becomes a product of factors which can not vanish in the usual range of 
the bound states \m\— m<i\ < Pq < (mi+m 2 ). The reduction condition becomes thus simply <I>(/j) ^ 0. The solutions 
$(po) °f the Bethe-Salpeter equation for which <I>(/i) =0 will thus not lead to solutions of the 3D equation. This 
would be the case for the solutions antisymmetrical in po, which may appear when the masses are equal. We have 
of course the possibility of performing a reduction based on a different constraint, like putting a fermion on the mass 
shell, as Gross |J . The Bethe-Salpeter amplitude is given in terms of the wave function by 

*(po) = V2^G (po)^r(po,M)*'- (27) 

At po = (J., we can factorize an operator (Pq — S)^ 1 out of Go(fi), the 3D potential V out of and use the 

3D equation ( |l9| ) to recover the relation (|25|). 



2.3 The QED case. 



To be definite, we shall work in the QED framework, but our results will be more or less directly adaptable to other 
interactions, and/or to the two-boson and boson-fermion systems Jl3| . The modified kernel Kt is then, in Feynman's 
gauge: 

ia 1 / ia \ f .. I I 



K T (p',p) 



47T 3 k 2 + i 

r 2ill5 



;(7i-72) + 



-Pn-5 



Air 3 J J 1 kl + iekl 
\ B Pip2 ) (71 -72) + 71*^01 {qi){ii-i2)G 02 {q'ic)l2n 



(71-72) Goi(gi)G 02 (<7i B ) 



(28) 



The various contributions are symbolized by the graphs of figure 1. In the reducible graphs, we use Gr — Gq — Gs 
instead of the Go we would use in the calculation of the scattering amplitude, so that the graph B, for example, 
symbolizes the sum of the second and third terms of (28). The external lines are there only to indicate the initial 
momenta pi,p 2 and the final momenta p'^p'2 of the fermions 1 and 2. The internal momenta of fermion 1 are qi,Q2, ■■■ 
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and the momenta of the photons leaving the fermion 1 line are k\, k 2 , .... All momenta will be written in terms of 
the fixed external momenta and of fci, ...k n -i, which we choose as integration variables. With these conventions, the 
different graphs of a given order differ only by the insertion points of the photon lines on the fermion 2 line. The 
internal momenta of fermion 2 will vary accordingly and we shall denote them by q'i,q' 2 , ■ with a supplementary 
index identifying the corresponding graph. The momenta of the intermediate fermion lines are thus: 

qi = pi - ki, q 2 = Pi - h - k 2 , ... all graphs. (29) 

q'i=P2 + ki, q' 2 =P2 + ki + k 2 , ••• ladder graphs. (30) 
nl — 1 permutations of ki,...k n in (30) crossed graphs. (31) 



3 One-body limit. 

3.1 One-body limit of the 3D reduction. 

We shall now take the one-body limit (7712 — >oo) of the 3D equation in QED. We shall first define E as a finite part 
ofP : 

P = m 2 + E. (32) 



For /j,, we have, in any part of (E 



m| E — TYl 2 



(33) 



2P 2 

This is also the asymptotic form to give to p' and po in ( p8| ) , because of the constraint. Combining this result with 
©.we get 

p'lO'Pw^E, p' 2Q ,p 2 o^m 2 . (34) 

For the intermediate energies, we shall define more practical variables. Considering, to be definite, the a 4 diagrams, 
we have 

qio=E-oJi, q 20 = E - uj 2 , .... (35) 

q'w = m 2 + u[ , q' 2a — m 2 + uj' 2 , .... (36) 

fcio = - 0, fego = wa - Wi, /c 3 o=w3-w 2 , fc 40 = 0-cj 3 , (37) 

ho = -(ho + k 20 + k 3a ). (38) 

In the ladder diagram (figure 1): 

w[ = kio, uj 2 = k w + k 20 , u' 3 = k w + k 2Q + k 3Q . (39) 

LU[ = LUl, LU 2 =0J 2 , UJ 3 — UJ 3 . (40) 

For the 4!-l other diagrams of order a 4 , we have simply to permute (kio, k 2 Q, k 3 o, k^o) in (|39|). This leads to an equation 
with only the /?2 = 1 components, where h 2 can be replaced by m 2 p3J: 



^ = 



1 



E — hi + iehi 



V'tp 



(41) 



V'(f,p) = -2in\^' + [i^) Idh 



la 
4^3 



47T 3 % 2 V 471-3 

dk 2 dk\ 



1 



1 



k 2 +ie kf + ie 



1 



LUl + if 



277T(5(cJl 



ll>i C + ie 



Goi(9i)/?i 



1 



1 



k 3 + ie k 2 + ie kf + ie 



1 



lu 2 + ie 



+ 2iir5{uj 2 ) 



1 



+ 277T(5(wi) 



with 



+--|G i(g 2 ),3iGoi(gi)/3i + 

Goi(<7i)/3i = -=; r / — * \ . ■ r 7=r: 

- - n-i(gi) + iehi{qi) 

qi = Pi - h, ?2 = Pi - ki - k 2 , ... 



(42) 

(43) 
(44) 
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Since lj' 1c = —lui, we see that the second-order contributions cancel mutually. It can be shown that similar cancellations 
occur at each order B [l0| O, |l2|, so that we remain with the first term of V' (a Coulomb potential) only. This 



result can also be obtained in other gauges. Equation (23) has then only one root, corresponding to the normal 



solution. At the ladder aproximation, however, we have only the terms in 

1 



\-2m8(ui) = (45) 

uji + it u>i — it 

with nothing to cancel them. Equation (|23[ ) can then have several roots. We have actually computed them at the 
zero-body limit ( m\ — >oo, static model) |lJj] , and we do not see how they could not be present for a finite m\. It must 
be noted that, if we make the instantaneous approximation without making the ladder approximation, the cancellation 
does not occur either: we remain with a dependence of the potential on the total energy, and possibly again with 
several roots of (p3|). 

3.2 One-body limit of the Bethe-Salpeter equation and its 3D reduction. 

In order to get the TO2^oo limit of the Bethe-Salpeter equation (Q), we shall write 

Po = m 2 + E, po = /J, - u. (46) 
Writing u>) = </>(w), K(p — Lu',/i — Lu) — k{uj' ,ui) we get, at the 7712—5-00 limit: 

= p z 1 T ■ h — ^~/3iM)H- (47) 
h — hi — lo + lehi lo + it 

The 3D reduction of ([47]) is obtained by writing 

= 1>{u>) + - - 1 . — di (48) 

h — n,i — u> + ieh\ uj — it 

= v^F*(w)^, ^ = F . 7 ftM(O). (49) 

Here, G# is obtained from Go by a simple change of the sign of one it (as it is also the case in Gross reduction [|J). 
The consequences of this change are not infinitesimal, as the integral of Gr with respect to lo is zero for the positive 
values of h\. Putting the iterations of in (p9h gives 



-2in 



E — hi + iehi 



U /.■(((.())+ I duiptK&ui)- 1 rvMwi.O) 

(— + ii — /ii + iehi)(LJi — it) 



+ I dw 2 duJi[3in{Q,uj2)- — 1 . — /3i k(w 2 , wi) 

[—0J2 + E — hi + ithi){u)2 — it) 



' / 9 1 «(wi,0) + 



(— wi + E — hi + ie/ii)(o;x — ie) 



V>' (50) 



With the exact kernel, the higher-order terms of k (crossed diagrams) cancel the higher-order terms of the expansion 
( pT)| ) (we get ( |4l| ) with the first term of fll^)). At the ladder approximation, the kernel is given by 

ia 1 

PiK{U2,q^;u}i,qi) = —» j-z — =; „ : . (51) 

47T 3 (w^-wiY - \qi — q\Y + it 



Equations (|5Q-|5lD give the same results as Equations (|41[-^2[) , both for the exact model and for the ladder approxima- 
tion. 

In order to perform this 3D reduction, we must have ip' — — iv / 27r(^0)(O) 7^0, (using ( f47j ) and (fi"9"|)). This implies 
that <j>(oS) must have a pole at origin (accordingly, ( ^5| ) shows that $(/x) must be infinite at the mi— >oo limit if we 
want to get a nonvanishing ). This is certainly true for the normal solutions, for which we can in principle compute 
tj/ by perturbations around a nonrelativistic approximation. We shall verify below that it is also true for the abnormal 
solutions of the ladder approximated Bethe-Salpeter equation at the zero-body (mi->oo) limit. 







4 Asymmetrical zero-body limit. 

4.1 Asymmetrical zero-body limit of the 3D reduction. 



The mi^oo limit of the reduction ( }4l}]42| ) will be performed after writing E = mi+W. At the ladder approximation, 
we get 

ia 1 



Wij)'{jf x ) = -2iir d A Pl 



47T 3 -(pi—fil) 



IC 



ia 

47T 3 



d 3 qi did 



1 ~2 



wf - (pi — qi) 2 +ie (— a»i + W + ie)(u>\ — ie) ujf - (qi-pi) 2 + it 



ia 

47T 3 



d qzdwzd > qxduji 



w 2 _ (P1-92) 2 + ie (-^2 + VP + ie)(w 2 - ie) 



1 



1 



1 



ip'(pi). 



(52) 



(u> 2 -wi) 2 - (<72-<7i) 2 + «e (— wi + W + ie)(wi — ie) ujf - (q\ -pi) 2 + ie 

With the exact kernel, we must keep only the first term. Since we expect to get the static model, we shall search for 
solutions in 

i>'{pi) = exp[ipi-f] (53) 

where x is a set of three continuus quantum numbers which can be interpretated as the conserved relative position 
of two infinitely heavy fermions. As the only remaining 3-momenta are k% = (pi — q\), &2 = (<fi — 9*2), etc.... in the 
propagators of the photons we shall perform the integrations with respect to the ki . Using the relations 



Pi = Pi~ (ki +k 2 + ...), 



3 exp[i k ■ x] —2ir 2 



d A k 



exp[i|cja;| 



(54) 
(55) 



removing a common factor exjp[ip\-x] and choosing |x| = 1, which means taking l/\x\ as the energy unit and 
rendering W and uj dimensionless (the energy will thus be W/|iS| from now on), we get 



W 



-2in 







—ia 


mi 


[2, + 





du>l exp[i|0 — uj\\ 



1 



2tt 



dio^dwi exp[i|0 — 
1 



{-uJi + W + ie){uji - ie) 
1 



— exp[i|cL>i — 0|] 



(— cJi + W + ie)(u)i — ie) 
With the exact kernel, we have simply W = —a. 



(-W2 + W + ie)(u} 2 - ie) 
— exp[z|wi — 0|] + . 



exp[z|w 2 -wi| 



(56) 



4.2 Zero-body limit of the Bethe-Salpeter equation and asymmetrical zero-D reduc- 
tion. 



The mi^oo limit of the Bethe-Salpeter equation (|47|) will be, after writing E = nii + W: 

1 



We shall search for solutions in 



(W - w' + ie)(uj' + ie) 

<f>(w,pi) = x(w) exp[ipi ■ x 



drpi duj k(u , Pl ;uj,pi)<j)(uj,pi). 



(57) 
(58) 



The only 3-momenta in k being again the ki of the photon propagators (in both the exact and the ladder approximated 
models), we shall again perform the integrations with respect to the fcj. Defining 



d 3 pi k(uj' , p 1 ^, iv , pi) exp(ipi ■ x) = exp(ip\ • x) 77(0/ ,uj). 



(59) 
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and choosing |x| = 1, we get 



At the ladder approximation, we have 

%Oi \ — %Oi 

k(u',p\;uj,Px) = -j-3 j-, r- 2 — 7^ -, , 2 , . , r)(u',u) = —exp(i\u'-w\). (61) 

Air 3 (to' — uj) z — {p\ — px) 1 le ^ w 

The zero-D reduction is obtained, as in section 3.2, by changing the sign of the second ie of (|60|), leaving a term in 
8{lo') as the initial term of an iteration. This gives ( fiff ) at the ladder approximation, and W = —a for the exact 
model. The condition for this reduction is 



duj 7?(0,w)x(w) ^ (62) 
which implies, via (|60|), that x(w) must have a pole at origin. 

4.3 Condition for the asymmetrical zero-D reduction of the ladder-approximated static 
model. 

We think that the absence of abnormal solutions of the 3D reduction of the exact Bethe-Salpeter equation at the 
one-body limit is probably the consequence of their absence for the original finite-masses Bethe-Salpeter equation 
itself. There remains however the possibility that such solutions exist and disappear as a consequence of the one-body 
limit, of the 3D reduction, or of the combination of both. It is of course impossible to study directly the behaviour of 
these solutions, the very existence of which is questioned. We know anyway that the abnormal solutions of the Wick- 



Cutkosky model survive the one-body limit 19 . We know also that the abnormal solutions of the ladder- approximated 
static model survive the zero-D reduction we performed in |l5| ]. This last zero-D reduction is however based on the 
choice u) = W/2 (just between the two poles at w = and lu = W) and treats the two fermions in a symmetrical way, 
while the present zero-D reduction is based on the choice to = , in order to get the very compact result W = —a in 
the exact model. We must thus check if the condition for this new zero-D reduction is also satisfied by the abnormal 
solutions of the ladder approximated static model. If we write 

W 

u(w)= X (w + y), (63) 



we get a more symmetrical Bethe-Salpeter equation 

1 

- W 2 - uj' 2 - ie 2ir 



= iw2 1 n ~ ~oT" / duj exp[i\uj'-uj\}u(uj). (64) 



4 

where the sign before the ie comes from the fact that we are interested in negative values of W (we had +ieW). 
Equations (^o|) or (|64|) are the Bethe-Salpeter equation for the static model, at the ladder approximation, in the 
relative energy representation. In a previous work |]l5| , we studied equation ( |64| ) in the relative time representation, 
performing a Wick rotation and a zero-D reduction. We could meet our previous results by performing a Fourier 
transform of ([34j), but it turns out to be also possible without leaving the relative energy representation. Writing 

U ^ = 1W2 1 ,2 ~ ( 65 ) 

jW 2 — uj 12 — ie 

W = -2w, lo = -wt, ^- = A, (66) 

ZTT 

we get 

... iX f , exx> \iw\t' — til , , ,„ , 

v(t') = - dt 2 ' " «(t • (67) 
w J t — 1 + ie 

The condition for the zero-D reduction ( p6[ ) is that v(t) must not vanish at < = — 1. If we derive twice equation (|6 
we get the differential equation 



v"(t) 



2A 2 
w 2 



t 2 - 1 



v(t) = 0. (68) 
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Although we are still in the relative energy representation, equation (|6^) and the complex conjugated equation (|67]) 
(written thus for v* ) are exactly the equations we studied previoulsly in the relative time representation |l5[ (the 
Fourier transform of the propagator takes the form of the potential, and vice- versa). We found that we had to solve 
the differential equation ([38]), and, using the Wick rotation in (p7|), that the solutions had to be square integrable 
along the imaginary axis. Equation (|6^) being regular and symmetric around the origin, the solutions are regular at 
origin and can be made even or odd. In fact, for each eigenvalue Wi, we found either an even or an odd solution, but 
not both simultaneously. A zero-D reduction of (p4|), with 

GoM = TWi 1 i ~ G s (w) = 5(u) f x ^ =-2nr^M ( 69 ) 

jW z — uj z — le J jW z — lo' 1 — le W 

demands u(0)^0, so that the odd solutions of the Bethe-Salpeter equation were not solutions of the reduced equation 
(we did not miss them, as they do not satisfy the normalisation condition deduced, like the homogeneous Bethe-Salpeter 
equation for the amplitudes, from the inhomogeneous Bethe-Salpeter equation for the propagator). 

To see if a solution of (|67]) leads to a solution of the reduced equation ( [j6| ) , we must verify that it does not vanish 
at t = — 1. The solutions v of (|68|), a Fuchs equation around this point, can be written as linear combinations of two 
solutions of the form 

oc oc 

« + = E C «( < + 1 )" +1 ' v- = Y J C-{t + l) n + cv+log(t+l) (70) 

n=0 n=0 

in the + < 2 domain of the t complex plane. At t = — 1, v + is zero, and v~ is Cq~, by definition a nonzero 
arbitrary constant. As written above, we know also that we must build even and odd solutions around t — 0. In order 
to build these solutions, we must in general use a linear combination of v + and v~ , which will differ from zero at 
t = —l, unless built with v + alone, which should then be even or odd when re-expanded into a series of powers of t. 
Although this would imply an infinity of cancellations between the coefficients , we did not exclude this possibility 
a priori. We checked numerically that it does not happen for the 20 first eigenvalues computed in ref. Jl5|| ■ 



5 Symmetrical zero-body limit. 

5.1 Symmetrical zero-body limit of the 3D reduction. 

We can take directly the m\ = m2 = rn— >oo limit of the 3D reduction. We have and keep fi = 0. We shall write 

P Q = 2m + W, (71) 

(72) 
(73) 
(74) 
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The fermion propagators, in the second-order terms of the exact model, will give a contribution proportional to 

1 



\W-wx 



(75) 



2" ' ' - 2 ' 

with uj' 1c — —u>i (at the ladder approximation, we would not get the crossed term) . In contrast with what happens with 
the asymmetrical zero-body limit, these terms do not combine to zero. In Coulomb's gauge, however, the only relative 
energy dependence lies in the fermion propagators, and the integral of ( f75| ) with respect to u>i is zero. The integral 
of the crossed term is indeed zero, as the poles lie both on the same (upper) half part of the wi-complex plane. The 
integrals of the two other contributions combine to zero (they correspond to the modified ladder term, Gr, which was 
designed for that). At higher orders, we find the products of the lower fermion propagators (W/2+ui' iG +ie)~ 1 , where 
uj' iG corresponds to the ith lower fermion propagator of the graph G. It can be shown that the ui' iG are combinations 
of the uii, with coefficients ±1,0. For the uncorrected ladder graphs, uj' iG —oJi- For the crossed graphs, at least one 
LOi appears only with minus signs in the lower fermion propagators, so that the integral with respect to this u>i will 
be zero. We remain then with the modified ladder graphs (built with the ladder term plus the 6 ) which give also zero 
by integration. We conclude that, in Coulomb's gauge, both the exact kernel and the ladder approximation give the 
normal solution only. This confirms our result above (and that of ||) about the absence of abnormal solutions in the 
exact static model, without contradiction with our result about the ladder approximation, which is not supposed to 
be gauge-invariant. 
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5.2 Zero-body limit of the Bethe-Salpeter equation and symmetrical zero-D reduction. 

Writing 



Po = 2m + W 7 po = —oj 



(76) 



in the Bethe-Salpeter equation at the ladder approximation, and taking the m — > oo limit, we get directly the 
symmetrical equation (64). The zero-D reduction of (p4J), using (K39), gives 



W = -2iix 







—ia 


(£)7 







1 



\W 2 -wf-ie 



-^j Jduiidw! exp[i|Q-«3|] | lw2 _^ 



9 , 

/ 



exp[i|w 2 -wi| 



2in 



exp — 0|] 



(77) 



The terms between brackets correspond to the two first terms (modified ladder part) of (75). 



6 Conclusions. 

We pointed out that the cancellation mechanism which transforms some well chosen 3D reductions of the exact 
Bethe-Salpeter equation into a simple Dirac equation (for two fermions in QED) at the one-body limit, implies that 
these reductions have no abnormal solution at this limit, nor at the zero-body limit (static model). Conversely, this 
cancellation does not occur at the ladder approximation, so that the abnormal solutions are present in this case. Our 
argumentation is based on the fact that the abnormal solutions, when present, are implicitly preserved by the 3D 
reduction, via the dependence of the potential on the total energy, and should thus disappear when this potential 
reduces to a simple compact potential, like a Coulomb potential. This cancellation mechanism in the 3D reductions 
at the one-body limit occurs in many situations (2 fermions, 2 bosons, 1 boson + 1 fermion with various interactions) 
|], |h], [ll], [l2|, [l3| Q , providing many new examples for the discussion of the abnormal solutions problem. The 3D 
reduction must be based on a constraint which puts, at the one-body limit, the heavy particle on its mass shell. 
Putting directly one particle on its mass shell, as in the 3D reduction of Gross without taking the one-body limit, is 
however not enough. These facts and the result of ref. || about the absence of abnormal solutions in the exact static 
model, support Wick's conjecture that the abnormal solutions are spurious consequences of the ladder approximation, 
and would disappear in the exact model. However, it remains possible that the exact model has abnormal solutions, 
which are killed by the one-body limit, the 3D reduction (if the corresponding Bethe-Salpeter amplitudes vanish for 
the relative energy fixed in the reduction), or the combination of both. This does clearly not happen with the normal 
solutions, and our study of the static model has shown that it does not happen either with the abnormal solutions at 
the ladder approximation. 

The main problem of the existence of abnormal solutions of the exact, finite masses Bethe-Salpeter equation remains 
thus open. What we have shown in this work is that the higher-order contributions to the 3D potential of well chosen 
reduced equations for the exact model cancel mutually at leading order in each of the masses. The dependence on the 
total energy, which is related to the presence of abnormal solutions, is thus much less important than at the ladder 
approximation, and possibly too small to allow for more than one root of Po-Ei(Po). It could perhaps be possible to 
eliminate this dependence completely order by order, using the large freedom which we still have in the choice of the 
3D reduction. We made clear anyway that, although the non-ladder contributions can be treated as correction terms 
in the study of the normal solutions, they modify completely (if not annihilate) the abnormal solutions obtained at 
the ladder approximation. The detection and computation of abnormal solutions goes thus beyond the validity of the 
ladder approximation. 
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